A natural origin of primordial density perturbations 



Richard Lieu^ and T.W.B. Kibble^ 



ABSTRACT 

We suggest here a mechanism for the seeding of the primordial density fluc- 
tuations. We point out that a process like reheating at the end of inflation will 
inevitably generate perturbations, even on superhorizon scales, by the local dif- 
fusion of energy. Provided that the reheating temperature is of order the GUT 
scale, the density contrast 6ji for spheres of radius R will be of order 10"^ at 
horizon entry, consistent with the values measured by WMAP. If this were a purely 
classical process, 5^ would fall as l/R^ beyond the horizon, and the resulting 
primordial density power spectrum would be P{k) oc with n = 4. However, 
as shown by Gabrielli et al, a quantum diffusion process can generate a power 
spectrum with any index in the range < n < 4, including values close to the 
observed n = 1 (5^ will then be cx for n < 1 and l/i?^ for n > 1). Thus, 

the two characteristic parameters that determine the appearance of present day 
structures could be natural consequences of this mechanism. These are in any 
case the minimum density variations that must have formed if the universe was 
rapidly heated to GUT temperatures by the decay of a 'false vacuum'. There 
is then no a priori necessity to postulate additional (and fine tuned) quantum 
fluctuations in the 'false vacuum', nor a pre- inflationary period. Given also the 
very stringent pre-conditions required to trigger a satisfactory period of inflation, 
altogether it seems at least as natural to assume that the universe began in a flat 
and homogeneously expanding phase. 



The large scale homogeneity and flatness of the early Universe as revealed by COBE and 
WMAP observations (Smoot et al 1992, Bennett et al 2003, Spergel et al 2007, Hinshaw et al 
2009) of the cosmic background radiation (GBR), together with the formation of the observed 
large-scale structures, present three of the greatest challenges to contemporary cosmology. 
The hypothesis of inflation (Guth 1981, Albrecht & Steinhardt 1982) offers an ingenious 
solution to these problems. 
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In this Letter we wish to focus our attention upon the supremely important problem of 
structures. An observational result of WMAP that may be interpreted as supportive of con- 
ventional models of inflation is the Harrison-Zel'dovich (HZ) power spectrum of primordial 
density modes in Fourier space, 

P{k) = 141' oc A;", (1) 

with n ^ 1. Provided P{k) 'levels off' at k ^ I/th where th is the comoving radius of the 
horizon, the density contrast Sr within a sphere of radius R is given by (see e.g. Gabrielli et 
al 2002) 

= ^ = / d^mikmk) ^^,n<l; and ^„ n > 1, (2) 

where u is the mean density of the universe, u + Su the mean density within the sphere, 
V — AnR^/S, and W{k) is the three dimensional Fourier transform of the top-hat window 
function W{r), defined as W{r) — 1 when r < R and W{r) — when r > R. 

The power spectral index of n f=:i 1, at least from the viewpoint of the ensuing 5\o^l/R^ 
scaling of real space density variations, is clearly indicative of a higher degree of spatial 
homogeneity than a Poisson distribution where n = and 5\ oc 1/R^. What is the origin of 
such an index? Zcl'dovich (1965), and Zel'dovich and Novikov (1983) were among the first 
to point out that small and causal displacements of particles from their original positions in 
a super-ordered state can only lead to an index of n > 4. Cosmologists then had to resort 
to more exotic explanations of n ~ 1, e.g. perturbations generated by quantum fluctuations 
during inflation. More recently, however, Gabrielli et al (2004) noted that the conclusion of 
Zel'dovich et al, while valid classically, is too restrictive in general: causally limited quantum 
interactions have access to the broader range, n > 0. This then opens the possibility of a 
novel and very natural kind of explanation for the observed spectral index n ~ 1. 

Let us first be clear about our choice of gauge here and after. Since superhorizon 
densities are ambiguously defined, eq. (2) is valid only in the synchronous gauge, where the 
amphtude evolution of the seed w.r.t. t occurs at the rate Sr oc r'fj{t) so long as R > rn 
and Sr <^ 1 (section 9.3.6, Kolb & Turner 1990). We are interested here in superhorizon 
evolution, i.e., for a given scale R, the period during which rn < R, between the epoch 
t = texit when it exited the horizon during the inflationary era, and the epoch t = tcntry at 
which it re-entered during the radiation (or matter) dominated era. As is well known, if the 
linear evolution equations are valid throughout this period, then the amplitude 6r will have 
the same value at tentry as it did at texit- At both limits Sr has physical significance and is 
unambiguously defined, although in the intervening period it is gauge-dependent. Moreover, 
this 'horizon-crossing' value of Sr also does not vary significantly with R. It is determined 
observationally by the low-harmonic temperature anisotropy of WMAP as Sr{R — th) ~ 3 x 
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10~^ (e.g. Bennett et al 2003), although the final four year COBE data favored the slightly 
lower number of 2 x 10~^ (Peacock 1999). Beyond the horizon, the observations tell us that 
the scahng of 5r is close to 



The constancy of 5ij(tentry) is conventionally taken to be a reflection of the very slowly 
changing conditions during inflation. 

In spite of the aforementioned milestone achievements, there remain signiflcant uncer- 
tainties and weakly justifled tenets in the inflation model, most of which have to do with the 
fine tuning of the initial conditions (see Padmanabhan 1993 for an account of many of them). 
Thus e.g. for every scale to exit with nearly the same amplitude, the Hubble parameter and 
the time derivative of the scalar field that drove infiation must be very precisely con- 
stant during the horizon exits of all density-contrast modes relevant to structure formation; 
indeed modifications of the basic scenario which lead to different patterns of spatial density 
fiuctuations do exist (e.g. Kawasaki et al 2003, Hall et al 2004, Yamaguchi & Yokoyama 
2004). Secondly, to trigger infiation the pre-infiationary universe must already be homoge- 
neous on supcrhorizon scales (Vachaspati & Trodden 2000). Thirdly, the formula for these 
delicately constructed zero-point fluctuations has both infra-red and ultra-violet divergences, 
and regularization is needed on each end. Fourthly, the result depends on a choice of initial 
state, usually assumed to be some kind of 'vacuum state', but this is somewhat arbitrary 
since there is no unambiguous deflnition of the vacuum in a curved background. Lastly, most 
scalar field theories of GUT phase transitions involve too strong a coupling constant; as a 
result, they yield an exit amplitude larger than 10~^ by orders of magnitude. Hence, the 
claim that infiation 'solved' the problem of structure formation by 'predicting' the observed 
density seeds may have been overstated. 

The purpose of this paper is to propose a way of improving the situation. The point to be 
demonstrated is that the origin of structures does not have to be the variance in the vacuum 
expectation value of at all: even for a completely classical, uniform the generation of a 
primordial density pattern with the observed amplitude is more than a mere possibility; it is 
perhaps even an inevitability. Central to our argument is a revisit of the 5i?(tentry) = ^^(texit) 
relation, which is normally seen as the reason why a classical vacuum cannot seed structures. 
That is, it is often argued that causal processes cannot create or modify density contrasts 
on supcrhorizon scales. Yet this is only approximately true. Particles generated during the 
reheating process can propagate some distance and so either penetrate or leave a chosen 
superhorizon-size sphere, though of course this can only happen for particles lying within a 
thin shell. Even if the energy density is perfectly smooth the reheating process will perturb 
all scales. 
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Nevertheless, if reheating were a purely classical process, then as Zel'dovich (1965) 
and Zel'dovich and Novikov (1983) showed, the rearrangement of energy could not yield 
a power spectral index less than 4. But reheating is not a classical process, and thus, 
according to Gabrielli et al (2004), is not so constrained, essentially because the correlations 
of quantum fields behave differently from those of classical fields. The real-space behavior 
of 6r described in eq. (2) or (3) is mathematically compatible with P{k) oc /c" for any n in 
the range 1 < n < 4. 

Now in the synchronous gauge the perturbed metric in comoving distance coordinates 
and conformal time is 



[dr = dt/a(t) with c = 1 for now). The Planck function for the energy density is, in this 
metric, the same as that in the unperturbed metric where hij = 0. To see this, we start with 
the mean and unperturbed occupation number n(p,T) = l/[exp{p^ /k^T) — 1] per quantum 
state of the primeval fireball, where p is the physical momentum, and = |p| the physical 
energy. When the metric is perturbed according to eq. (4) , we may use the results of Ma & 
Bertschinger (1995), who found that the momenta conjugate to are Pi — a{Sij + hij/2)p> 
(and also Pq = ~flP°)) a-nd that the mean of the momentum-energy density tensor is 



where Tq = aT, (-g)-^ = a-\l - h/2) and dPidP2dP^ = a^{l + h/2){pydp'^dn, with 
h = ha, where dQ is the clement of solid angle. It is evident that to lowest order the effect 
of the metric fluctuation li cancels, so that e.g. in the perturbed metric. 



which is the same as that in the unperturbed background Friedmann space-time. 

From eq. (6) ensues a 'benchmark', or reference, level of spatial fluctuations for the early 
radiation universe, viz. the spatial thermal noise of an equilibrium photon gas of volume 
V and radius R (both comoving), with standard deviation 6^^ = (16A;b/3S')"'^/^, where S 
is the entropy in this volume. Since the entropy is conserved during the later adiabatic 
expansion of the universe, this is also the entropy in the same volume today, namely S = 
(43/22) X (AuoV/Tq), where Uq is the energy density of a photon gas at Tq — 2.725 K, and 
the factor of 43/22 arises because the entropy today includes not only that of photons but 
also that of neutrinos. (The entropy of the photons today is equal to that of photons plus 
electrons and positrons before pair annihilation. Hence the ratio of the photon entropy to 



ds'^ = a^(r)[-(ir^ + {Stj + h.ij)dx''dx^] 
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neutrino entropy now is (2 + 4 x 7/8)/(6 x 7/8) = 22/21. See for example Peacock 1999, 
pp. 279-281.) Thus we find 



Note that because this only depends on the entropy, which is conserved, eq. (7) does not 
involve the extra degrees of freedom that may have arisen as a result of the proliferation 
of particle species in the early universe when temperatures were very high. Thus in the 
context of thermal noise the gas at any such epoch may roughly be treated as an ensemble of 
randomly moving photons, with Planckian energy spectrum. This then affords the following 
heuristic derivation of eq. (7). Observe that the mean photon energy e — hV and its standard 
deviation 5e are both ~ /cb^o, so that for a volume of total photon number N and total energy 
E ~ Nk^To the variance in E is {5Ef = N{5ef ~ N{kBTof, or 

(SEV _ /6uy _ 1 _ 1 / he y / f'^ x^dx 

\¥) ~ \^o) ~ ~ 8^ [k^oj [Jo e--l) ' ^ ' 

where the integral in parentheses equals 2.4. The square root of this is only 15 % larger than 
the correct answer of eq. (7). 

It is important to note that the part played in structure formation by the perturbations 
of eq. (7) has not been sufficiently recognized. This is most easily seen by applying the 
hnear growth equation (eq. 5.123, Peebles 1993) to a small subhorizon region i? <^ r^j, for 
which expansion is negligible and the gravitational term may be ignored. In that case it 
reduces to the equation of undamped sound wave propagation, predicting zero growth. We 
know, however, that on such small scales, where local thermodynamic equilibrium (LTE) 
can quickly be restored by particle diffusion, any super-thermal fluctuations are damped, 
and any sub-thermal ones are enhanced, in relatively short times. 

During the reheating phase between the cosmic times U and tf a, classical inflationary 

vacuum would have dissipated into heat. Since the energy of this 'false vacuum' is uniform 
at ti, the density contrast in superhorizon volumes will still be small at tf. Nonetheless, the 
important point is that the contrast cannot be exactly zero because thermal energy is not 
vacuum energy. Once particles are present, then on the surface of each such volume there 
exists a thin layer within which particles located outside this volume can random walk into 
it, and vice versa. The comoving thickness of this layer, assuming 6t = tf — ti is not too 
large, so that we may treat H{t) as nearly constant and equal to its radiation-era value at 
tf, is 

Z"*^ edt eSz c{zi — Zf) 2 
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where g is the number of helicity states of light particles at the time of reheating (fermions 
counting with a factor 7/8), and fl^ is the normalized radiation density parameter (including 
only photons, not neutrinos). Thus the energy variance is [SEY = N'{5eY ~ N'{kTY where 
A^' 3N6R/R. The effect, as already explained, is neglected by the linear perturbation 
theory. Since the total energy remains at ^ ~ NkT we now have 



where the last step was taken with the help of eq. (7). 

In this way, density seeds with the scaling relation of eq. (2) became available quite 
naturally. Specifically, 5^ oc l/R^ corresponds mathematically to n = 1 + a where a 
is arbitrarily positive. Physically, as we remarked earlier, small diffusive movements of 
particles in the quantum-interaction regime can yield a power spectrum with P{k) oc /c" and 
< n < 4. In the case of n < 1, eq. (2) says that Sj^ oc l/i?^"*"" for R > r^/, the scaling with 
R now differs from that of eq. (10). This is not surprising, because eq. (10) was derived 
assuming a classical horizon, which does not always lead to the same i?-dependence of 6^ as 
the scenario of the quantum particle diffusion. The recent WMAP data (Hinshaw et al 2009) 
give as the best fit n the slightly lower value of n = 0.96, or Sj^ oc 1/R^-^^, very close to the 
form of eq. (10). 

Turning next to the question of amplitude, the factor that determines the significance 
of the role of diffusion. As we have emphasized, reheating is a quantum process, so ideally 
a full quantum treatment would be required. However, while a classical calculation may not 
necessarily reveal the correlations inherent in a quantum process, they should give a reason- 
able estimate of the magnitude of the effect. In reheating, the universe made a transition 
from vacuum energy domination to radiation domination, during which rn was changing 
slowly (from a decreasing function of time to an increasing one). The value of r^ for this 
period is determined by the reheating temperature Treheat, which we assume to be of order 
the GUT scale, i.e. one could adopt 



Tr\ C I 2 / A B^ielieat \ 



rH~ ^ = —^ _ _r^pa 200 cm. (11) 

Oreheat-f^reheat ^reheat Hq y gflj \10 ^ GbV J 

We may then use eqs. (9) and (10) to obtain 

Sn - 1.29 X 10-^ ( "-^) V- V for R > rn. (12) 
^ VlO^GeVy VlOO/ \ ) yrn) ~ 

It is reasonable to assume that {zi—Zf)/zf < 1 (which just means reheating is not protracted) , 
but that this ratio is not very small compared to unity. In that case, we see that 5r for the 
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horizon entry scale ol R — th is the amplitude of fa 10~^ observed by WMAP and required by 
structure formation models, i.e. eq. (2). By comparing eq. (3) with eq. (12), one sees that 
both the desired shape and normalization of Sr can emerge from this mechanism. 

For each superhorizon region enclosed by a real sphere the surface diffusion will continue 
with time beyond t — tf as 5r oc t^^'^ during the radiation era, since 5^ oc 5R oc (eqs. 
(9) and (10) with a{t) oc t^l"^). This occurs independently of the gravitational growth of 
linear theory, viz. 5r oc rfj cc a'^(t) oc t in the synchronous gauge, which commences ait = tf 
(when Tfj fully exhibits its radiation dominated time dependence) and persists until re-entry. 
Since the latter scales with time much faster than the former, there is no need to take into 
account the 5r oc t^/^ rate at epochs t > tf. The role played by the diffusion process is really 
to 'plant the seeds' during the period ti < t < tf when linear growth had yet to get under 
way, i.e. if initial density contrasts were absent by t — tf, there would be no seeds available 
for developing into the structures of today. 

In conclusion, the mechanism we proposed will certainly generate primordial density 
perturbations, and under very reasonable assumptions these could have the required ampli- 
tude and spectral index to agree with WMAP observations. There is then no need to postulate 
additional quantum fluctuations in the vacuum to solve the 'problem of seeds'. As we have 
stressed, this is possible only because reheating is an inherently quantum process. At the 
present time we do not have a specific quantum model to propose that would have the desired 
properties. Rather we assumed that such a model exists, and pursued the consequences of 
this assumption. If the proposal stands up, then the development of a specific working model 
will be an important task for the future. If the ideas presented here can actually account for 
all the essential parameters by the time the data have the quality to clinch them, the search 
for a simple and viable inflationary hypothesis might be freed from constraints and greatly 
simplified. 

The authors thank Ruth Durrer, Francesco Sylos Labini, and Robert Crittenden for 
helpful discussions and the provision of vital technical information. 
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